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A MATHEMATICAL MUNCHAUSEN 


Paul Schillo 


At the time of this writing, ten years have elapsed since the un- 
timely death of 1. Neustadt. During this decade, the reform movement 
1a methematics, which he pioneered, has been marking time - almost as 
if it were afraid to proceed without a leader of his imposing stature. 
But, in order to imbue mathematics with Neustadti an methodology, what 
is needed most is not a leader with Neustadt’s revolutionary zeal but 
a few more revolutionary discoveries like he was able to make. 

A sufficient reason for the failure to discover new mathematical 
truths in the manner outlined by Neustadt is the lack of disciples in 
his school. Every reformation must expect the forces of reaction to 
have their innings, but, in this case, the game is being forfeited 
without contest. Most mathematicians have an erroneous conception of 
Neustadt’s position, some have had no occasion to use his famous law, 
with the result that they give him little thought, and instances are 
known where otherwise mature mathematical scholars have failed to 
acquaint themselves with even the rudiments of his philosophy. There- 
fore, while this article is simply a brief account of the man and his 
major work, it is hoped that it might serve to correct some of the 
misconceptions of Neustadt’s intent and so help to proselyte his 


program. 
One notion has it that Neustadt sought to destroy mathematics. 


Nothing could be further from the truth; his orientatia was positive 
at all times, for what he wanted to see was its reconstruction along 
modern lines. His sole aim was to make the structure of mathematics 
functional, to have it adopt scientific method and, thereafter, take 
its place among the modern sciences. Because of the cruel way in which 
hasty critics have distorted his progressivism, it is a wonder that 
he has not been accused of having suggested the tearing down of the ca- 
thedral of Notre Dame before the erection of the United Nations’ building. 

A good example of the way in which he has been misunderstood is 
the popular interpretaion of his facetious disproof of Malthus’ Law. 
From the readily verifiable fact that two can live as cheaply as one, he 
proceeded to show by conventional mathematical methods that ten billion 
can live as cheaply as one. To naive critics, this might seem to be an 
attempt to discredit mathematics as a whole; actually, it was intended 
merely to show the absurdity of mathematical induction, a technique of 
which he strongly disapproved. He had no quarrel with fundamental, 
time-tested devices like addition and multiplication and used them 
frequently in his experiments. 
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It was not without justification that Neustadt came to the con- 
clusion that mathematics was suffocating in its steadily worsening 
state of introversion, He felt that it had accepted the egocentric 
predicament centuries ago and had given up trying to extricate itself. 
To him, it was autopophagy, all this waste of cognitive energy in pro- 
ducing theorems, which were no more than reworded statements taken 
from the particular axioms that happened to be in vogue. “The present 
pace of mathematics is that of a snail,” he once said. “And what is a 
snail? It’s just a snake that has swallowed its tail” It was this 
irritation with all self-centered philosophical positions that goaded 
him on to his classic proof of the existence of a third personal re- 
ality,. the essence of which is found in his often-quoted remark: “Er 
iszt; daher ist er!” He was sure that mathematics had been feeding 
off its own fat for so long that it had become easy prey to the pre- 
datory natural sciences which picked at its bones like it was cold 
turkey. He envisioned its rebirth and hoped to teach it how to hunt 
for new truths in order to survive in a fiercely competitive age. 

Unfortunately, the complete story of Isack Neustadt has yet to be 
written. The information about his early life is fragmentary. If he 
wrote anything prior to the formulation of the great law which bears 
his name, there is no record of it. Only a sketchy description of his 
formative years has been gleaned from the close-mouthed villagers of 
Dilldorf, Pennsylvania where he was born and raised. His buddies in 
the Grand Amy of the Republic, in which he served as a private until 
he suffered a head wound at the battle of Gettysburg, have added Little 
to the portrait, except to agree that he kept to himself and was af- 
fectionately referred to as “Pop”. 

Concerning his life in Dilldorf after the war, this much is known: 
The little indications of budding genius caused the villagers to con- 
clude that he had been hexed, a consequence of having taken up arms 
and of having been encamped for a long time near the wicked city of 
Philadelphia. He began to read non-religious books, some of them in 
English and some of them concerned with the heinous heresy of science. 
Also, he showed no inclination toward physical labor and did nothing 
but pick a little fruit during the harvest in order to get money for 
books. It is to the lasting credit of the gentle people of Milldorf 
that, despite their ignorance of I sack Neustadt’s greatness, they saw 
to it that he did not want for anything essential to his survival. He 
was given shelter in the smithy and shared the victuals in the biack- 
smith’s dinner pail. In fact, he Lived more comfortably than any of 
his benefactors. But the need that these good people failed to satisfy 
was his greatest need, the hunger to be understood. The admirable 
patience, with which he endured those pangs throughout his early years, 
was to be given good use later when he had to face the hard fact that 
most mathematicians were simply Milldorfers with academic degrees. 

His escape from the subtle tyranny of his idyllic, native village 
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was one of those remarkable little accidents which, taken all together, 

play the dominant role in the pageant of science. During the roaring 
twenties three gangsters out of Cicero, Illinois rubbed out a Wet 
lobbyist in Washington. The federal men soon knew who the hoodlums 
were 48 well as their devious getaway route via the backroads of five 
states. But, until they spoke to Isack Neustadt, none of the witnesses, 
who had seen the black limousine on its chicken-killing, wash- dusting 
journey through the hinterland, could help them refute its ower’ s 
claim that it had not been outside of the Chicago environs since its 
purchase. And it was only by a stroke of luck that they were able to 
contact Neustadt. For, due to their dislike of the government and its 
minions, the people of Dilldorf treated the investigators coldly; and 
those who spoke English maintained with righteous indignation that, 
during the fruit harvest none of them could possibly have been so 
idle as to gape at a “devil's chariot”. But, just as the government 
men were about to drive off, a village wag made some remark about “der 
verdammte Faulpelz”, which caused the dour Dilldorfers to laugh most 
heartily. Questioned about this outburst, an unusually helpful young 
interpreter explained that they were amused by the thought of one who 
might not have been so busily engaged at the time of the vehicle's 
Passing, and that was I sack Neustadt. 

As it happened, Neustadt had a vivid recollection of the incident. 
In keeping with his neighbors’ opinion of him, he had been napping 
under an apple tree alongside the road. He told of having a strange 
dream about the war, a part of his life about which he had all but 
forgotten. A falling apple struck him on his temporal scar and he 
awoke in the same state of intense alertness as when he had regained 
consciousness in that gully at Gettysburg. It was at that moment that 
he saw the car pass by in a swirl of dust. He could not make out any 
of its details, not even its color, but its license plate seemed to be 
suspended motionlessly before him. It was an Illinois plate with the 
number 2-]2-21-2. Even if he had not been endowed with a photographic 
memory for numbers, he would have remembered that particular one, 
because he had a profound respect for numerical palindromes. 

Since his testimony was critical in the successful prosecution 
of the three torpedoes, he received a five thousand dollar reward 
which, after taxes, was reduced to $2,122. 2 by a remarkable coincidence. 
However, by going to court and by accepting money for just remembering 
a number without doing any laborious reckoning, he became an anathema 
to the people of Dilldorf. Therefore, he hiked to Harrisburg and hopped 
a freight train on itsempty run to Florida. He took up residence in 
a boathouse on the ocean side of “liami Beach, where, because of his 
frugality, he was able to live out his life comfortably on the re- 
ward money. It is a matter of public record that he never received 
remuneration for his many published works; money meant nothing to him. 

In the stuffy refrigerator car which took him south, he began to 
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contemplate his great law. Meditating on the events which had altered 
his life, he considered the cause and affect relationship that had 
to exist between the falling apple and the remarkable number 2- 12-2 1-2, 
He thought of another falling apple in the past and of the man who 
had seen its significance. After all, was not the latter’s name merely 
the anglicization of his own? A combination of his native mysticism 
and the semi-delirium, induced by the stifling atmosphere of the car, 
led him inevitably to the conclusion that he was the reincarnation of 
Isaac Newton and that there was something of great importance to be 
discovered in the number 2-]2-21-2. Whether or not one is disposed to 
believe that there was anything supernatural about such an accumulation 
of coincidences, it must be agreed that the subsequent behavior of 
Isack Neustadt constitutes an excellent illustration of the James - Lange 
theory. For, even if he were not Isaac Newton, he acted enough like 
his prototype to make any distinction seem irrelevant. 

Hoth literally and figuratively, it can be said that Neustadt’s 
Law was incubated in the warm sand of the Atlantic Ocean, for, during 
his first year of peripatetic study on the beach, he wrote the number 
1-1 2-21-2 over and over again on that golden slate. When the first 
burricane in his experience seemed about to rip the boat-house from 
its pilings, his most brilliant idea struck with greater intensity 
than the wind. He converted the dashes in 2-12-21-2 to minus signs 
and put plus signs before all but the leftmost of the other digits, 
obtaining 2- 1 +2- 2+ 1 - 2, a finite series with zero sum. He 
knew that other six-term series had the property of adding to zero, 
but he was not to be dissuaded from the further consideration of this 
particular one. 

Within a month from that fateful moment, he made the final altera- 
tion in his series, whence it became .2 - 1 + 2- 2+ 1 - .2. which 
as the reader knows, is now called“ Neustadt’s basic series” or,by 
some writers, simply “Neustadt’s base’. According to his own account, 
the first arithmetic progression which he used was 1,2,3,4,5,6. When 
he produced the derived series (or derivative) .2- 2+6-8 +5- 1.2 
by multiplying the successive terms in the base by the respective 
terms in the progression and observed that it too had a zero sum, he 
was so overwhelmed by emotion that he had to sit down beside a palm 
tree and there fell into a state of unconsciousness. 

However, with the discovery that, for any six-term arithmetic 
progression, the derivative of his base had zero sum, his work had 
just begun. He was then faced with the herculean task of proving that 
his observation did indeed have the status of a mathematical law. 
This he accomplished in seventeen months without any outside assistance. 

Before beginning the experiment which was eventually to prove his 
law, he took as his control group twenty zero sum bases different 
from .2- 1+ 2-2 +1- .2 and listed two thousand distinct arith- 
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metic progressions to be used on the bases. Among these progressions 
there were many like 2.71828, 3.14159, 3.56490, 3.98821, 4.41152, 
4. 83483 so it should be evident that he did not limit the generality 
of his results by selecting special progressions which would make 
his work easier. His ability to calculate mentally helped him to 
complete this experiment in a shorter period of time than it would 
have taken ten ordinary mathematicians. But, despite this demonstra- 
tion that an exceptional person can undertake an immense mathematical 
project and, single-handedly, bring it to a successful conclusion, it 
ought to be clear that the search for mathematical truths can best 
be carried on in well-equipped laboratories by highly trained staffs. 
In the future machines and teamwork must be at the heart of mathematics 
or it shall forever remain an agglomerate of results from uncoor- 
dinated individual effort. No doubt, those years spent in solitary 
labor on his famous law led Neustadt to organize the first mathematical 
cooperative and to preach the doctrine of centrally controlled research 
as an alternative to the inherently feeble laissez faire mathematics. 
But his hopes that this small group would be the nucleus of a uni- 
versal organization were based on the faulty assumption that others 
were as devoted to this single objective as himself. The cooperative 
lasted less than four months and did nothing constructive with its 
modest resources except to contribute to the Lemke-for-president 
campaign fund, 

On the first day of his experiment to show that, for any six-term 
arithmetic progression, the derivative of .2-1+2-2+ 1- .2is 
itself zero sum, Neustadt had a fortunate accident. Using the arith- 
metic progression -3,-1,1,3,5,7 on his base, he had just obtained the 
zero sum derivative -.6 + 1+2 -6+5- L4, his 45th verification 
of the law without a single contradiction. He was getting tired and 
mistook this derivative for one of the bases in his control group. 
When he applied the progression -3,-1,1,3,5,7 to -.6+1+2-6+5- 14 
he found that its derivative 1.8 - 1 + 2 - 18 + 25 - 9.8 was zero 
sum, something that had not happened previously to a derivative of a 
base in his control group. And, the next morning, when he realized 
his mistake, he discovered something more amazing. By applying the 
progression 1.6,1,.4,-.2,-.8,-1.4 to 18 -1+2 - B+ 25 - 9.8, a 
derivative of a derivative (or a second derivative) of his base, he 
obtained a third derivative with zero sum: 2.88-1+.8+3.6- D+ 13.72. 
He immediately revised his program for the experiment in order to 
prove the law as itis formulated today: For any arithmetic progressions, 
all derivatives of any order of the base .2-1+2-2+ 1- .2 are 
zero sum. 

Right after this law was copyrighted on Flag Day in 1932, many 
mathematicians thought that it could be proved much more easily by 
their conventional devices than by Neustadt’s scientific method, How 
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wrong they were is evidenced by the fact that, despite numerous attempts, 
none of them have been able to come up with a deductive proof of it. 
Their sole contribution has been a forceful demonstration of the fact 
that the time Neustadt devoted to his successful project is negligible 
when compared with the time used in their futile efforts to belittle 
his methodology. 

Now, unlike many mathematicians, Neustadt was no absolutist.He 
knew that the real numbers were phenomena like atoms and organic cells 
and that the laws which govern their behavior can be ascertained only 
to a degree. This degree of certainty, or truth value, of a known law 
must be quantified by assigning to it a number between zero and unity, 
commonly called its probability. He determined the probability of the 
truth of his law to be a number which differed from unity well beyond 
the thousandth decimal place, which is so close to certainty that, 
should it ever fail, one might do well to look for some catastrophic 
alteration in the universe. The manner in which he obtained this pro- 
bability from his experimental data depended on a special application 
of the famous formula: 


a + 
Zz = 
n 
which Euler communicated to Diderot. As if anticipating the current 
popularity of electronic computers, he used the binary system because 
of its simplicity and, consequently, set a = -2 and 6 = 2, He let v 
be the number of experimental verifications of his law and c the 
number of contradictions. Then, n was made | or the vth prime according 
as v was or was not zero. After he evaluated the integer x, the re- 
quired probability was determined by the empirical formula: 


P = 2°¢ Q°2e°8 


Because his experiment yielded no contradictions of the law, this 
formula was reduced to the simpler formula: 


P=1-2% 


and P approached 1 asymptotically as the experiment progressed and 
the number of verifications increased. 

With the completion of this work, Neustadt did not rest on his 
laurels. But it is not the intention of this article to go into all 
that he did. It is well known that, as a champion of science, he led 
many crusades. Yet, the power of leadership failed to corrupt him and 
he retained throughout his life that golden streak of humility which 
is the banner of true genius. He was as cordial with the mountain folk 
of central Florida as with the polished aristocrats of suburban New 
York. So unobtrusive was he in his frequent contacts with people along 
the beach that, unless they had gone to Miami with the intention of 
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meeting him, they were unlikely to realize how important a person he 
was. This aura of simplicity, which made him such an amiable man, in- 
fluenced even his most profound thoughts. It enabled him to see that 
complex numbers were just ordinary numbers perceived by complicated 
people. 

It is vain but interesting to contemplate what new truths I sack 
Neustadt might have uncovered had he not met with sudden death in 
his favorite retreat by that seaside coconut palm, under which 
he would sit by the hour like the great Gautama. Who knows what secret 
of the universe was locked up in that last message which he left in 
the sand? Perhaps someday a devoted disciple of his will look at that 
eleven term mathematical legacy: 


.2-4+ 27 - 96+ 210 - W2.4+ BA- B2+ 81 - 20+ 2.2 


and will see its significance. But, as a sagacious humorist put it at 
a recent testimonial dinner, it should be hoped that this favored 
fellow will be lying under an apple tree and not a coconut palm when 
he gets his inspiration; otherwise, he might not live to develop it. 


Buffalo, New York. 
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Electrical Analogies 


a new field for mathematicians and physicists 


In applying electrical analogies to various physical phenomena, 
mathematicians and physicists at Lockheed have developed a highly 
effective method of solving complex dynamic and static problems. 


They have been establishing the analogies through use of a Network 
Analyzer so advanced in concept that it is one of, perhaps, three or four 
instruments of its type in the country. Their activities involve solving 
partial differential equations by substituting finite difference equations. 


Through these techniques, they have been able to solve problems 
that are, at best, extremely difficult to solve on even the largest 
and most modern digital computers. 


Most of Lockheed’s endeavor in this field lies in unexplored areas. 

The dynamic and static problems are on a scale never before approached. 
Both creativity and ability to exercise individual initiative are required 
because of the frontier nature of the work in evolving the analogies and 
relating them to a wide range of phenomena. 


It is a new and ideal field for those possessing a physics major and 
mathematics minor or mathematics major and physics minor. 


The growing importance of this area has created a number of new 
positions. Because this type of work is so recent in origin, experience is 
not required. Inquiries are invited. Address E. W. Des Lauriers 

at Department 90-15, 


Advanced Education 
California Institute of Technology Lo k h ec ec d 

It is expected that those joining Aircraft t orporation 
Lockheed will have opportunities to 
pursue this subject in Lockheed-paid C alifornia Division 
courses at California Institute of Burbank 


Technology. The first course is being 


offered in the spring semester. ( A al fo rn t a 
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THE ISOPERIMETRIC PROBLEM IN THE PLANE 


Glen E. Bredon 


1, We shall consider here the isoperimetric problem for convex 
regions in the plane and its extension to more general sets. Part I 
is completely elementary as regards the methods used and presupposes 
next to nothing on the part of the reader. In it we prove the existence 
and uniqueness of the solution to the isoperimetric problem in the 
plane, restricted to convex regions. The uniqueness of a solution, if 
one exists, can quite easily be extended to regions bounded by a Jordan 
curve, but the existence becomes a little more difficult to show. In 
the second part we consider the existence proof for these more general 
regions, and also formulate and prove a very general isoperimetric 
theorem valid for arbitrary measurable sets in the plane. This part 
presupposes a very little knowledge of point set topology and measure 
theory, and a few statements are left unproved. Some of the places 
where the words “clearly” or “obviously” are used are admittedly actual 
gaps in the exposition. The author feels, however, that these gaps can 
be easily filled in by the intelligent reader by straightforward, if 
perhaps long, reasoning. 


Part I. 


2. A convex region means here a closed set & in the plane such that 
if p and q are in R, then the line segment joining them lies in R. 

The 1soperimetric theorem for convex regions in the plane says that 
among all convex regions in the plane of a given area A, the circular 
disc of area A, and only this disc, has minimum circumference. Clearly, 
this is the same as saying that the circular disc has the largest 
area among convex regions with fixed circumference, and again, these 


are equivalent to saying that for any convex region of area A and cir- 
C 


cumference C, A S C*/4n, with equality if and only if the region is 
a circular disc. (A and C are understood to be finite. ) 

We note here for completeness that the length of the boundary of 
a convex region, and indeed, in a certain sense, the length of any 
curve, is defined as the least upper bound of the lengths of the in- 
scribed polygons. Note that if C is finite then the lengths of all 
the chords, and hence of those from any particular point, are bounded 
and so the convex region is itself bounded 
Lemma 1. If y andy’ are convex curves (i.e. curves that enclose convex 


. 


regions) with lengths C and C', and if y surrounds y’ then C2 C 


Proof. By the definition of C' we may assume that y' is a convex polygon. 
The proof in this case is effected by successively replacing arcs of 
y by line segments which are just the sides of the polygon y' produced 
until they meet y (see fig. 1). This operation clearly decreases cir- 
cumference and we finally arrive in the end at the circumference of 
y', giving us our result. 
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Figure 


This lenma shows, among other things, that every bounded convex 
curve lias a finite length. (That is, it is “rectifiable’.) We need 
only enclose the curve in a suitable circle. 


Next we develop briefly the interesting operation of symmetrization, 
due to Steiner. (See (5) and also Pélya (4).) 


3. Given a bounded convex region K and a line l, the “axis of symmetri- 
zation”, we define K, the symmetrized image of k with respect to l as 
follows. Consider lines r which are perpendicular to | and cut #, Let 
r' be the segment of r which is the intersection of r with R. Let r, be 
the segment of r which has the same length as r’ and as its midpoint 

on Ll, We then define Hh, as the union of all these line segments r ;. 
(See Fig. 2). It is easy to see that Rh, is again convex. 


Figure 2 
Lemma 2, Area (,) = Area (R) and Circ (R,) < Cire (A). 


Proof. That the area is invariant is an elementary application of 
integral calculus. We reduce the second statement to the corresponding 
one for polygons, Let P; be a polygon symmetric to / and inscribed in 
R, such that Cire (P,) 2 Cire (K,) - €, where € is anarbitrary positive 
number, Let P be the polygon inscribed in HK whose vertices are mapped 
by the symmetrization into the vertices of P,. P, is then, in fact, 
the symmetrized image of P with respect to l. If the theorem is true 
for P and P; then we have, using Lemma 1, Cire (K,) < Cire (P)) + € S$ 
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Cire (P) + ©« < Cire @) + © and, since € is arbitrary, the theorem 
would be proved. The theorem for polygons clearly reduces to showing 
that the sum of the lengths of the lateral sides of a trapezoid with 
fixed bases and altitude is least when the trapezoid is isosceles, 
which in turn reduces to showing that among triangles with common base 
and given altitude the isosceles one has the least perimeter. This 
latter can be seen by use of an ellipse whose foci are at the ends 
of the base of the triangles and whose minor semi-axis is equal to 
the common altitude of the triangles. Since the sun of the distances 
from any point on an ellipse to the foci is constant we have (see 
fig. 3) AC + CB 2 AD + DB = AE + EB giving our result. 


D> C 


Figure 3 


Lemma 3. I1fR is a bounded convex region and not a circular disc then 
there exists a centrally symmetric region /}' which is convex, not a 


circular disc, of the same area as Kk and of no greater circumference 
than 


Proof. We note first that any figure F, all of whose synmetrized images 
are circles, is itself a circle. If suffices to show that, for any n, a 
regular 2” — gon can be inscribed in F. Select an arbitrary axis of 
symmetrization l and let P, be the chord of F which goes into the 
diameter of Ff, perpendicular to |. This is our 2 - yon. Say that a 
_ gon has been inscribed in Then, given a side S of 
consider the two symmetrizations in lines respectively perpendicular 
and parallel to S. The first shows there is a point of the boundary 
of F at the correct distance from S, and the second shows that this 
point must be precisely in the correct position to be a vertex of a 
regular 2" — gon builton the vertices of P60 This discussion is 
admittedly vague and the reader is urged to make a drawing and con- 
vince himself of the validity of the argument. 

We now return to the region R in the statement of the present lemma, 
By the above we can assume that A is symmetric with respect to some 
axis |, since we can always find an axis of symmetrization which doesn’t 
send & into a circle. Now diviae # into two parts by this axis of 
symmetry and flip one side over as in figure 4. This clearly turns R 
into a centrally symmetric region with the same area and circumference, 
and which is convex and not a circular disc, 
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4. Proof of the [soperimetric Theorem. 


Theorem 1, (Uniqueness) If among all convex regions of circumference 
C, h has a maximum area, then A is a circular disc. 


Proof. Using lemma 3, we know that if A is not a circle then there 
must be a convex, centrally symmetric region R’ with the same area and 
circumference as Hh which is also not a circie. Let g be the greatest 
radius of R' and s the smallest. Then by using the construction in- 
dicated in Figure 5 we see that s = x, for, if not, we have obtained 
a region with area equal to (y - s)* sin y + area (R) and the same 
circumference as RK, But the only centrally symmetric region with all 
radii equal is the circular disc, which shows that R must be a circular 
dise after all. 


Figure 4 Figure 5 


Theorem 2, (Existence) The circular disc has a maximum area among all 
convex regions with given circumference. 


Proof. Say that K is a convex region with circumference C and area A 
such that C*/A < 4m. We can assume R to be centrally symmetric by 
lemma 3 and hence can apply the operation indicated in Figure 5, which 
clearly gives a region which is again convex and centrally symmetric. 
Let HK, be the region which results from n applications of this operation 
and denote by g, and s, the greatest and smallest radii of K, respec- 
tively and by Y, the acute angle between them. Each application of 
the transformation leaves the circumference C fixed and increases the 
area by (yg, - s,)? sin Y,. Since RK, 1s convex and s, the smallest 
radius, K, must lie completely to one side of a perpendicular to s, at 
its endpoint. Hence we see that s, 2 4, cos¥, and thus 


3 
(g, - s,)? siny, -s,) /6, 


Since the area cannot increase indefinitely (it is bounded, for in- 
stance, by (C/2) 2) we must have (4, ~ s,) ~ 0. Now consider two circles 
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of radius s, and 4, inscribedin and circumscribed about R, respectively. 
Recalling lemma I we have the inequalities 


2nms, S2mg, and 78) SA, 


which show that s, and ¢, approach a quantity r with the property 


that C = 2nr and Lim A, = A' =n. But since A,2 A we have that 


4n = G*/A'S C?/A < 47, a contradiction, proving the theorem. 


Part Il. 


5. We have proved the isoperimetric theorem for convex regions. The 
uniqueness part of the theorem is easily seen to hold for regions 
bounded by a Jordan curve since any such region which is not convex 
clearly cannot be a solution to the problem. We will indicate two 
methods of extending the existence proof to these regions. Hoth methods 
are concerned with the so-called convex hull ACA )of a revion R. This 
is defined as the set of all points lying in triangles with vertices 
in Rk. When K is closed, (fh) is convex and is clearly the smallest 
convex set containing Rh. Obviously the existence of a solution to the 
isoperimetric problem for the regions we are considering will follow 
if we can show that the circumference of A(X) is not larger that that 
of K. (It is easy to see that the area is actually greater unless R&R 
is already convex.) 

The first method of showing that cire (h(h)) cire lies in 
observing that the boundary of A(R) consists of subares of the boundary 
B of & together with straight line segments joining points of B and 
also that these straight line segments replace arcs of B joining the 
same points. Anyone fairly familiar with point set topology should 
not have much trouble proving this statement and the existence proof 
of the isoperimetric problem follows immediately. 

The second method rests on an interesting expression for the length 
of an arbitrary rectifiable curve, which we state as a theorem but 
will not prove, although the proof isnot difficult (see Blaschka (1)). 


Theorem 3. For any rectifiable curve y with length L we have 


@ 
(*) L = ¥ n(p,v) dp dy (Lebes que) 


0 


where n(o,v) is the number of intersections (possiblyinfinite) of y 
with a line of polar angle ¥ and polar radius op. 

Granting this we have immediately that circ (h(f)) < cire (HK) since 
the corresponding n(?,¥%) satisfy this inequality except possibly at 
the set of measure zero wiere the first of these is infinite. 
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6. By use of formula (*) we can state the isoperimetric problem for 
the plane in a very general form as follows. Given an arbitrary measur- 
able set M in the plane, we define the circumference of M by the 
integral (*) where n(o,/) is now the number of intersections of the 
line (p,¥%) with the boundary of © (‘= complementation, 

= closure). 


Theorem 4, 1f M is a measurable set in the plane with area (Lebesque 
measure) A and circumference C, then A < C?/4n, Furthermore, if 
M= int (W) (i.e. if M is the closure of its own interior) and C is 
finite, then equality holds only if Y is a circular disc. 


Proof. Since the inequality is trivial in case C is infinite, we will 
assume C to be finite. Given M we will make a series of replacements 
which do not decrease the isoperimetric quotient A/C? and will finally 
arrive at a convex set by these replacements, from which the first 
part of the theorem will follow by previous results. 


i) Area (M2 Area (M and WNW CH NM’ go that Circ < Circ 
We let R = W, 


ii) Let A = int &),. (Note that now K = int (K).) This operation 
is easily seen to be the result of dropping a certain set of boundary 
points, so that the circumference is certainly not increased. Say 
that the area is decreased. This implies that there is a set B of 
boundary points for which m(B) > 0 (Lebesque measure). Now for fixed 
we have mB) = de, where m(p,¥) is the linear Lebesque 
measure of the set of points in which B intersects the line (¢,¥). 
Since this is positive, the set of » for which B intersects the line 
(o,/) in an infinite number of points is of positive measure. Since 
this is true for all ¥, we have that C is infinite, contrary to assump - 
tion. Hence the area of K equals that of R. 


iii) We define a relation ~ on int (K) by: p ~ q (p,q in int (K)) 
if and only if there is no k in K such that p and q lie in different 
components of kK - {k}. This is clearly an equivilence relation on 
int (K) and since the equivilence classes are open sets there can be 
at most a countable number of them, Let {E,} be the set of these 
classes and set F = and F F . Note that K F consists of 
boundary points alone and so by the argument under ii) does not con- 
tribute to the area, Now each F, is met by the rest in at most a 
countable number of points and hence Area (K) = Area (F) = = Area (F,) 
and Cire (kK) 2 Cire (F) = > Cire (F,). Let 


a = lub {Area (F )/(Cire (F )) 2} 
n 


(whose existence will follow later). Then we have 
Area (K) = Area (F)= = Area (F) <a> (Cire (Fy)? < 


a Cire (F,))? =a (Cire (F))?< a@ (Cire (K))? 


| 
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with strict inequality in the case that there are more than just one 
F_. Hence there is an F, such that 


Area (F) Area 


(Cire (F))* (Cire (F,))? 


iv) Now consider a single KF. We see that if any line meets F in 
exactly one point, then F, must lie completely to one side of the 
line, and hence by formula (*) the circumference of the convex hull 
n(F,) is less than or equal to that of F.. We have thus arrived at a 
convex set whose lsoperimetric quotient is not less than that of ¥, 
and hence, by part I, we have proved the first part of our theorem. 
Previous remarks and the essence of the above proof render the second 
part of the the theorem obvious. 


7. I take this opportunity to point out to the interested reader 
some of the different and highly elegant treatments of isoperimetric 
problems existing in the literature. Especially recommended is the 
excellent treatment by Polya (3) and also that of Gourant and Robbins 
(2). “ost of the groundwork was done by Steiner (5), who, however, did 
not bother much with the question of the solutions existence. Some 
more advanced material is to be found in Blaschke (1), while Polya (4) 
gives some highly interesting applications of symmetrization to mathe- 
matical physics. 
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The Second in « Series of Announce 
ments on Progresswe Expansion of 
Program and Mathematics 
at the Knolls Atomic Power Laboratory 


The steadily advancing nuclear program 
at Knolls Atomic Power Laboratory calls 
for new and imaginative departures in 
mathematics — ranging from the most ab- 
struse formulations of fundamental prob- 
lems to the digital solution of physical 
problems. To meet the consequent expan- 
sion of its Mathematical Analysis Program 
the Laboratory plans to increase signif 
cantly the number of qualified mathema- 
ticians now at work here 
openings have been created, in fact, to 
more than double the present mathema- 
tical staff. Mathematicians at all degree 
levels are invited to join this expanding 
program 


enough new 


As previously announced, a modern build 
ing is now under construction, principally 
for the use of mathematicians and phys 
cists. This Center will be equipped with the 
finest of facilities, including digital com- 
puters that rank among the most powerful 
available” Here mathematicians, working 
both independently and in association with 
theoretical and expermnental physicists, will 


WILL RECEIVE IMMEDIATE ATTENTION 


SCHENECTADY 


GENERAL ELECTRIC'S 


KNOLLS ATOMIC POWER LABORATORY 


MODERN MATHEMATICAL CENTER 


A LETTER TO DOR & ACKER EXPRESSING YOUR INTEREST 


Knolls Alomie Power Laboralory 


OPERATED FOR AEC BY 


GENERAL @® ELECTRIC 


IT 1S NOW DOUBLING THE 
STAFF OF MATHEMATICIANS FOR ITS 


enjoy an atmosphere in which the creative 
mind may find its full fruition 


As members of the Mathematical Analysis 
Unit, they will participate in the formula- 
tion of theories to describe new physical 
uituations now bemg encountered, in eval- 
uating these theories and adapting them 
to numerical solution by digital computers, 
and in evaluating reactor designs. Design 
evaluations will focus on the calculated 
behavior of mathematical models and will 
employ the most modern techniques in 
computer programming. The nature and 
complexity of these operations call for cre- 
atively new approaches and fundamental 
advances in these techniques. These mathe- 
maticians will also have the opportunity to 
deal with basic research in physics, chem- 
istry, metallurgy and many other aspects 
of nuclear science 


The program at the Knolls offers the atmos- 
phere, the equipment, the richness of sub- 
ject matter and the material benefits con- , 


ducive to a satisfying career in applications 


of mathematics 
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STABILITY OF PERODIC 
TIME-VARYING SYSTEMS 


Louis A. Pipes 


Summar y 


This paper presents a mathematical analysis of the stability of a 
class of mechanical systems that contain time-varying parameters 
of a periodic type. The mathematical analysis leads to the solution 
of linear second-order diflerential equations of the Mathieu-Hill 
type. For purposes of illustration the general methods presented are 
applied to the analysis of the motion of the “inverted pendulum’ . 


I. introduction 


The general theory of oscillatory mechanical systems and their 
electrical analogues whose parameters are linear and constant has 
been extensively developed and is well understood. In recent years, 
considerable attention and effort has been directed to the analysis 
and to studies of the performance of dynamical systems whoue para- 
meters vary with the time. Special cases of technical importance in 
mechanical engineering are described in the treatises on mechanical! 
vibrations by Timoshenko (1] and Den Hartog [2]. Time-varying 
systems of importance in electrical Engineering and communications 
include the induction generator, frequency-modulation circuits, the 
condenser microphone, and the super-regenerator. The theory of these 
circuits is summarized in the treatise by Mclachlan (3). 

The behavior of this entire class of dynamical systems can be de- 
scribed by differential equations that by suitable transformations 
can be reduced to an equation of the following general form, 

d*x 
(1.1) — =- F(¢iz = 0 
at? 


where the function F(t) 1s a single-valued periodic function of 
the time t of fundamental period T. 

This second-order linear differentia! equation is known in the 
mathematical literature as a differential equation of the Mathieu- 
Hill type, (3). The genera! solution of (1.1) depends on the nature 
of the periodic function F(t). If this function can be represented 
by a Fourier series of the form, 


@ 
F(t) + A cos (nwt ) + Bsin(nwt) 


n=l nel 
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where w= 27/T, then (1.1) is known as Hill’s equation. If, on the 
other hand, F(t) can be represented by the first two terms of the 
Fourier series (1.2) so that it takes the form, 


().13) Fit) = Ao + A, cos(wt ) 


then (1.1) 1s known as Mathieu's di%erential equation. 


The classical! solution of Mathieu's diflerential equation leads 
to linear combinations of Mathieu functions (3]. The classical 
solution of Hill’s equation involves computations with infinite 
determinants [3]. In general, the classical solutions are difficult 
to apply to special problems of practical importance. 

It is the purpose of this paper to describe two methods for the 
solution of equations of the type (1.1). Une method involves matrix 
miltiplication and the second method is an approximate one originally 
given by lord Rayleigh, [4], but which has been used very little 
in the technical literature of vibrations. Although these methods 
are general, in order to provide a concrete physical example of 
their use, they will be applied to the study of the motion and 
stability of the “inverted pendulum, (2,5). 


II. Equations of Motion of the Inverted Pendulum. 


Consider the inverted pendulum 
of Figure }. This pendulum con- 
sists of a weightless rod of 
length L to which is attached a 
mass m. The base of the pendulum 
is moved with a vertical accelera- 
tion S by a crank or cam mechanism. 
It will be assumed that the pen- 
dulum undergoes a retarding 
frictional force at its pivot point 
0 of « viscous friction type. If 
the angle @ that the pendulum makes 
with the vertical 1s smali, it can 
be easily shown 2 that the eq - 
uation of motion of the pendulum is 


. 1 
(2.2) +2b6- L (g+S)@ = 0 


The term 260 represents the retar- 
ding effect of the viscous frictional 
force of the prvot and g is’ the 
acceleration of gravity. In order Figure 1. 
to remove the first derivative term 


Inverted Pendulum 


a 
\ 
/ 
| 
/ 
| / 
| 
6) 
S(t) 


1956) STABILITY OF PERIODIC TIME VARYING SYSTEMS 


introduce the following transformation, 


d* 
dt L 
This differential equation is of the Mathieu-Hill type (1.1) with, 


Two cases of the motion of the base of the pendulum will now 
be considered. In the first case, it will be assumed that by means 
of certain cam mechanism the base of the pendulum undergoes the 
following square-wave acceleration, 


T/2 


(2.5) = K, 
< T 


S =-K,- 1/2 


< 
< 


where K ys a constant. 


In the second case it will be assumed that through the action 
ofacrank mechanism, the base of the pendulum is given the following 
harmonic motion. 


(2.4) S = a cos(pt) 


The effect of these two distinct types of motion of the base of 
the pendulum on its behavior will now be analyzed in detail. 


III. Motion of the Pendulum Produced by a Square-Wave 
Acceleration of its Base. 


In order to analyze this case, let, 


(3.1) 1 =p + = - K) 


The diferential equations of the motion of the pendulum (2.3) 
now take the form, 


(3.2) x - = 0, 17/2 
(3.3) = 0, 


Let the initial angular displacement and angular velocity of 
the pendulum at t = 0 be OG, and wo respectively. These initial 
conditions in © lead to the following initial valuesiin x(t) by 
means of the transformation (2.2). 


(3.4) x(0) = =x, (0) 
(3.5) x(0) = x,(0) 
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In order to determine the subsequent motion of the pendulum 
starting with the initial conditions (3.4) and (3.5) it is necessary 
to cascade the solutions of di™erential equations of the form (3.2) 
and (3.3) over intervals in time of duration T/2. This can be done 
by the use of a matrix mltiplication scheme [4]. 

In order to do this, the solution of (3.2) may be written in the 
following matrix form, 


x,(t) cosh (k,t) sinh (k ,t)/k, x ,(0) 
(3.4) = ° 

kysinh (k,t) cosh x (0) 
with the notation x, = x and Xo = x. 


Let 


(3.7) 6, = k,T/2 


then at the time t 17/2, the solution of (3.2) may be written in 
the form, 


x (T/2) k, sinh (6) cosh (6,) x ,(0) 

The solution of (3.3) may also be written in a matrix form 
similar to (3.6) with k, replacing ky. The initial conditions of 
the solution of (3.3) are the values x, (T/2) and x,(T/2) given by 
(3.8). Hence the solutions for x, and x, at t = T may be written 
in the convenient matrix form, 

(3.9) 


x (7) kS(6,) C(6,) k,S(6,) x,(0) 


(3.10) =k, C(6) = cosh S(6@) = sinh (6) 


If the indicated matrix multiplication is performed, (3.9%) may 
be written in the following form, 


x, (7) 


x,(T) 


— 
| A B x ,(0) 
(3.11) 
ic x, (0) 
where, 
- 
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cosh (6,) cosh (6,) + (k /*,) sinh (6,) sinh (6,) 
sinh (6,) cosh (6,)/k, + sinh (6,) cosh (6,)/k, 
k (6,) sinh (6,) + k sinh (6,) cosh 


= cosh (6,) cosh + (k sinh sinh 


.The solution at the time t = nT or after n complete cycles of 
the variation of the acceleration of the base, can be obtained by 
casading the solution (3.11) n times, It is, 


x, (nT) A Bi” x, (0) 


(3.13) 
x, (nT) x, (0) 


As a consequence of the transofrmation (2.2), the angle of the 
pendulum at the time t =nT is given by, 


(3.14) O(nT) = x (nT) 


The stability of the pendulum therefore depends on the behavior 
of x, (nT) for large values of n. By means of Sylvester's theorem 
in matrix algebra [46], it is possible to write, 

A Bi" (S.,,7DS,) (BS, ) 
(3.15) 

C 1 (CS,) 7AS,2 
where 


(3.14) S.= sinh fan) 


n 
(3.17) 

cosh(a) = 5 = cosh(6,) cosh 6,) + + k5/2k sinh 6, )sinh 6,) 


Therefore as a consequence of (3.12) and (3.14) the angle of 
the pendulum at ¢ = nT may be written in the following form, 


baT 
(3,18) 6 (nT) = [x - DS,) + *,(O)BS,) 


Therefore it can be seen that for large values of n the angle 
© will increase indefinitely and the pendulum will become unstable 
in the vertical position if, 


(A+D) 


(3.19) a > bT, where a = cosh"! 
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If, on the other hand, 
(3,20) a < bT 


the angle of the Pendulum 6 will eventually approach zero and the 
pendulum will be stable in the vertical position. In the absence of 
damping, 6 = 0 and the vertical position will be an unstable one 
unless a is a pure imaginary number in which case the hyperbolic 
functions in (3.18) become trigonometric functions. In order for 
this to be possible, it is necessary that, 


+D 
(3.21) 


If (3.21) is satisfied, the pendulum will perform oscil lations 
about the vertical position. Since ky is always positive, it is 
necessary that k, assume negative values in order for the condition 
(3.21) to be fulfilled. This means that the acceleration K must 
exceed that of gravity. 


IV. Harmonic Motion of the Base of the Pendulum, 


A study of the motion and stability of the inverted pendulum 
when its base is forced to execute a harmonic motion of the form, 


(4.1) S =a cos (pt) 


by the action of a crank mechanism will be undertaken in this sec- 
tion. The method of analysis is an approximate one suggested by 
Lord Rayleigh, [4]. 

If (4.1) is substituted into (2.1), it is seen that the equation 
of motion of the pendulum in this case takes the following form, 


(4.2) 6 + 2b6 - ; - ap = 0 


To simplify the form of this equation, introduce the following 
change in variables, 


(4.3) A(t) = 'ty(t) 
and 
(4.4) x = pt 


In the new variables x and y, (4.2) takes the following form, 


d’y 
+ b?)/p? -a cos (x) = 0 


(4.5) 
Now let, dx 
(4.4) w* = (y/L + b?)/p 


(4.7) =afL 


4 

) 
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Then (4.5) may be written in the convenient form, 


(4.8) 2 cos (x)] y 


This is a form of Mathieu's di%erential equation and its general 
solution has the form 3, 


=+@ 
(4.9) y= (nok )x 


where y = (-1)* and the a, quantities are constants. & is a para- 
meter to be determined by requiring (4.9) to be the solution of 
(4.8). If (4.9) is substituted into (4.8) and like powers of e/™ 
are equated, the following set of homogeneous equations in the a's 
are obtained, 


(4.10) [(n + k)* w*)a, - cla, +a ) 


-1 n+l 


where n runs through all values from - infinity to + infinity, and 
use has been made of the result, 


(4.11) 2 cos (x) = e/* + 


In order for the homogenous set of equations (4.10) to have 
other than the trivial solution a, = 0, for all values of n, it is 
necessary for the determinant of the coeficients of (4.19) to vanish. 
The vanishing of this infinite determinant determines the value of 
the parameter k. The recursion formula (4.10) may then be used to 
determine the a’s in terms of two arbitrary constants. This is an 
outline of the classical Floquet theory of the solution of Mathieu's 
equation {3}. 

If, however, the parameters of the pendulum are such that 


(4.12) 
(4.13) 


as is usually the case, then the following approximate method for 
the solution of (4.8) given byLord Rayleigh [4] my be used. 


Instead of the infinite series (4.9) the following three term 
approximation for y is assumed, 


neel 
(4.14) = a ei 


This approximation is justified because as a consequence of (4,12) 
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and (4.13) the values of the a’s of higher indices are very smal} 
and can be neglected. As a consequence of the smallness of w* and c? 
it can also be shown [4] that the parameter k is also of the order, 


(4.15) ke? <<] 


Hence the recursion formula (4.19) gives the three equations, 


+ cay = 0 


(4.14) ca_, tw*)a, + ca,= O 


Ca 4a, = 0 


where use has been made of (4.12) and (4.15). In order for this 
system of homogencous equations to have a non-trivial solution, the 
determinant of the coeficients of the @’s must vanish. Hence, 


-] c 0 
(4.17) c -(k?4w?) 
0 c 


If this determinant is expanded, the following equation for the 
parameter k is obtained 


(4.18) k? = 2c? y? 


If 2c? > w?, then (4.18) gives two real roots for the parameter 
k of the form, 


(4.19) k =+ (2c? ~ w?)? 


Fquation (4.14) then gives the following two linearly indepen- 
dent solutions of (4.8), 
Ye = +1+ce/*] 
(4,20) 
= +1 + ce/*] 


where C, and C, are arbitrary constants. Since the equation (4.8) 
is linear, the general solution is the sum of these solutions and 
it can be written in the form, 


(4.21) y= ** + [ce"J* +1 + cel*] 
This can be written in the alternative form, 
(4,22) y =A cos (kx + p) [1] + 2 cos (x)) 


where A and @ are two new arbitrary constants. In the absence of 


’ 
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damping, b = 0 and the angle of the pendulum may be obtained from 
(4.3) and (4.22) in the form, 


(4,23) =A cos (wot +p) [1+ cos (pt) 


where 
(2.24) 


The condition of stability is that the parameter k be a real 
number and hence 


(4,25) 2c* =w 
or 


(4,24) p > 


That is, in order for the pendulum to execute stable oscil lations 
of the form (4.23) about the vertical, the angular frequency of the 
oscillations of its base must satisfy the inequality (4.24), 

If damping is included, then the pendulum will perform damped 
oscillations about the vertical position of equilibrium of the form, 


(4,27) 6(t) = Ae” cos (wot +) [1 + cos (pt)] 


provided, 
(4,28) p> (2g, + L*b?)* 


where 


(4,29) » =|—>- 
L 


V. Conclusion. 


The two methods presented in this paper for the analysis of 
dynamical systems whose equations of motion are of the Mathieu-till 
type are quite general and useful in cases where the classical 
theory is dificult to apply. The matrix multiplication method is 
particularly useful when the variable parameters are of the square- 
wave or saw-tooth type. The approximate Rayleigh method is a 
very useful one when it can be applied and seems to have been neg- 
lected in the technical literature. 
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AN ELEMENTARY APPROACH TO THE USE OF THE 
RATE OF CHANGE CONCEPT FOR SOLVING PROBLEMS 


Raphael T. Coffman 


Foreword 


This paper presents a method for deriving expressions for the rate of 
change of certain algebraic and trigonometric functions and for applying 
the rate of change concept to derive area and volume formulas without, 
at any stage of the development, introducing the use of limits. 

From the time of the ancient Greeks all recognized methods of deriving 
such formulas, except for the simplest, require the use of limits in 
some form, ranging from the Greek method of exhaustion to the integral 
calculus. 

The writer has not found, in searching through such writing on 
mathematics as were available to him, any indication that this method, 
or any other not involving the use of limits, has previously been used 
for solving the type of problems considered in this paper. Some state- 
ments were encountered to the effect that for specific problems no such 
method is known. 


1. Introduction 


A line which moves in a direction perpendicular to itself may be 
considered to generate an area by its motion, the rate of productionof 
area is the product of the length of the line and its velocity, Fig. I. 


Fig. 1 Fig. 2 


If the line of motion is not perpendicular nor parallel to the line 
the velocity vector will have normal and parallel components, vx and 
vy Fig. 2. It is apparent that moving a line in a direction parallel to 
itself does not produce area, but only changes the position of the 
line. In Fig. 2 the area swept out by tie motion of | in the direction 
of vector v in any time will be a parallelogram, the area of which is 
1 x altitude. The altitude will be proportional to vz, regardless of 
the value of v, for equal time of motion. Hence: 


1. The rate of area generation by the motion of a line ts equal to the 
length of the line and its normal (perpendicular) velocity. If the 
line is changing in length, one end will have a parallel velocity. 
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This velocity does not generate area, hence cannot affect the rate at 
which area is being generated at any specific time. It determines the 
rate at which the rate of area production is changing, assuming uniform 
normal velocity. 

In what follows, all rates of change will be given in terms of one 
uniform velocity, which may be defined, without reference to limits, as 
change per unit time. If a velocity varies and its magnitude can be 
related to some other quantity, the expression showing the relationship 
sufficiently defines tne velocity. Thus if vy = 4xux, vy varies with 
x and vx and is specified as fully as is A in A =7r? for any value 
of r regardless of whether r is a constant or a variable. 

In the problem to follow the proposition (1) given above and the 
following are the basis for most of the proofs: 


2. If two quantities, originally zero, vary in such a manner that their 
respective rates of change are always in a constant ratio their 
magnitudes will always be in this ratio. 

This is true for any uniform rates of change, hence must be true 
if the rates of change vary, for it is true for whatever rate of change 
the quantities may have at any instant and they can change only to 
another value for which it is also true. 

In the notation used the rate of change of a quantity will be in- 
dicated-by v, as vx, v@, etc; where coordinates are used vx will be 
assumed to be a uniform rate of change and vy wll be expressed in 
terms of ux e.g. vy = 3x vx, 

The area bounded by a curve, the x axis and the ordinate at x will 
be designated as Ax, similarly with respect to the y axis for Ay. 

Volumes may be considered as generated by moving planes and the same 
arguments apply as given for generation of areas by lines. Hence: 


3. The rate of volume production by a moving plane is the product of the 
area of the plane and its velocity perpendicular to the plane. 


II. Problems 


A. Consider the area Ax in the first 4 vy v 
quadrant, Fig. 3, under the line 

y = x. By geometry this area is known 

to be 4x*. This may be found by the 

rate of change method as follows: Ay 


vAx = yux, which since y xX 1S XUX Ax WA 


vAy = xvy, vy = vx, hence vAy = xvx x 


[he total rate of area production = 
is xux + xux = 2xux and the area Figure 3 
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produced is xy = x? hence vx?=2xv0x 


Ax 
Ax+Ay  vAx+vAy 2Qxvux 2 


Ax = 4(Ax + Ay) = 4x? 


Rh. Next consider the area Ax in the 
first quadrant under the curve 


2» 
y = x*.The total area A is xy * 
xx? = vy 


vAx = yux = x?ux AY 


vAy = xvy, vy = vx? = Qxux 


from problem (A), hence, 44 | Vx 
Ax 


=A 


Figure 4 


fhe formula has been derived for the area under the parabola y = x” 
and an expression obtained for the rate of change of x’, namely 3x‘ux. 
by following the same procedure as for y = x’ the area under the curve 


y= x? is found to be % x” and the rate of change of x" found to be 


4x?vx. Thus it can be shown that the rate of change and the area under 
any curve of the form y = x" are respectively nx"*! vx and 
[1/(n + 1)] x"*! where n is any integer. 

The same procedure may be used to 
show that the same formula applies 
when nis a fraction of the form 
l/a where a is an integer, 


C. The equation y = /2, Figure 5 


is the inverse of y= x‘, hence = 

Ax (2/3)¥?. The same result may 

be obtained as follows: The rate Ay 

of change of x* has been shown to be ' 

2xuux; x may be regarded as (/x) ?, Ax Vx 

ux v(/x)? vx from which 


Figure 5 
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vy= fu ux 


vA = vAx + vAy = ux +% vx = (3/2) Vx vx A = =x'-5 


Ay  vAx Je tx 2 
vA (342) Vx vx 3 


Ax = (2/3)A = (2/3) x**° 


Both v(/x) and the area under the curve are of the forms 


- x"! 


n+l 


%y following the same procedure for y =x'/? y = x*?°, and so on it 


can be shown that the same general equation is true for any fraction of 
the form n = l/a_ (a = integer). 

To find the rate of change of a product, such as xy, construct a 
rectangle whose sides are x and y, Figure 6. 


N 
Wy 
Figure 6 Figure 7 


Then vA = yux + xvy = vxy. If A is constant yux + xvy = 9, since a con- 
stant has zero rate of change. 


1). Find the formula for the volume of a pyramid, cone or any pyramidal 
prism, in terms of the area of the base and altitude. 
The rate of change of a cube is vx? 3x*vx (Pig. 7). The side xis 
increasing at a uniform rate. 


If the side of a square in- 
creases at the same uniform rate 
as the cube in Fig. 7 and has a 
velocity vx perpendicular to its 
surface, the rate of volume pro- 
duction is x’vx and the solid 
produced is a pyramid with base x? 
and altitude x. [By giving the 
generating square the desired com- Figure 8 
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ponent of velocity parallel to its surface, the shape of the pyramid 
may be made regular or oblique as desired. A rotary motion in the plane 
of the square xx may be imparted, producing a spiral pyramid, 

The rate of change of the pyramid, x’vx, is 1/3 that of the cube, 
hence when x of the cube and pyramid are equal the volume of the pyramid 
is 1/3 that of the cube or (1/3)x?. vx 


Figure 9 Figure 10 


In Fig. 6 an area, A, is shown in the generating plane xx. This area 
has a constant ratio to the area of the square, that is A/x’ = constant. 


vA A ux 


Hence ——>= constant = ->—. The volumes produced will be in the same 
vx x ux 


Volume of solid A 4 


ratio. —— = -> from which, Volume of solid A 
x 


3 
3 A*x where A is the areaof the base and x is the altitude of the solid. 


E. Find the formula for the rate of change of the sine, cosine, tangent, 


etc. 


1. To find the rate of change of the sine and cosine construct, as shown 
in Fig. 10, a circle with radius 1, and tangential vector vx which re- 
presents the rate of change of the angle x. AB = sin « and O6 = cos x. 
The vertical component of vx is cos x vx and, since it is the velocity 
with which A is moving upward, it is the rate of change of AB. Hence 
v sin x = cos xux, Simlarly, the vector sin xvx = the rate of change of 
cos x and since cos x is decreasing this change is negative. 
v cos x — Sin xX UX, 


2. To find the rate of change of the tangent and secant, draw Figure II. 
Here since r 1,AB tan x and (A = sec x, 
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= gec xX UZ 


The vector AC is the rate of change of 
the tangent since this is the velocity 
with which the length of AB is changing. 

AC 
— = sec x or AC = sec x v, 
v 
1 
but v, ~ sec x ux 


AC = sec x+ sec x vx = sec*xux 


2 


v tan x = sec*xux 


v, is the rate of change of sec x (OA) 


Vo =v sec x = tan x + sec x Ux = 


sec x + tan xvx, 


3. To find the rate of change of the co- 
secant and cotangent: 


In Fig. BP, 


= CSC X UX 


rate of change of OA 


cot x 
Figure 12 


= cot xz = cotx sec xX Ux 


Since v, acts in a direction to decrease OA, the rate of change 
is negative, 


v csc x = — sec x cot x ux 
AC is the rate of change of AB and is negative because AB is decreasing. 
AC 


— = cse x or AC = csc x vy 
v 
1 


‘A 

v y 

= tan x or v, = tan x a) 

| = 
| Figure 11 
OA of 
‘2 
Vx 
V> 


v cot x == csc x vy = - esc x csc x vx = - esc’ x vx 


1. Find the rate of change of sin” x. 

In Fig. 13, make OD = AB = sin x, then DG = sin*x. 

The point D has two components of velocity a normal, DE = cos x ux 
since the rate of change of sin x = cos xvx and a tangential, Dil = 
sin xvx, since the tangential velocity is proportional to the distance 
from the center, 0. The rate of change of sin*x is the vertical com- 
ponent of the resultant of the two velocities. This component is 


KJ = KE + EJ 

KE = DE sinx =cos x sin x vx 
EJ = DH cos x =cos * sin X ux 
KJ = 2 sin x cos x ux 


This result can be obtained in a less complex fashion by considering 
a square whose sides are sin x. The rate of change of the area, which 
is sin’x is 2 sin x cos x vx. 


Figure 13 Figure 14 


F. In the preceding derivations of volume and area formlas the ratio 
Ux 
am has been constant, with the result that the desired quantities 
were obtainable in terms of this ratio, When the quantities are not in 


a constant ratio to each other a modified procedure is necessary. 


In Fig. 14 it is desired to find the area under the curve y = f(x) 
between the values of x = 0 and x = a, 

This area is equal to the area of some rectangle whose base is a 
and whose altitude is d, that is Ax = ad. d then is the average height 
of y in the interval under consideration. (Notation: y = average y.) 
The average value of vy is given by ~ a vy = - va. This does 

vx a a 
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not directly assist in finding the required area though it does suggest 
a method of approach: if another equation, y, = f,(x), for which 
vy, = f(x)(vx can be found, the average value of vy, will be the 
average height of y. Let f,(x) = ¢ when x =a and f,(x) = h when x= 0, 


vy, =~ y 


a 


- f = Ax 


As an illustration of the method, find the area under any portion 
of one arch of the sine curve from x = 0 to x = a. 
The rate of change of cos x is — sin x vx hence f, (2 = - cos x. 


A = cos a (-cos 0) = 1- esa 


If the equation of the curve contains more than one term the area 
under the curve will be the sum of the areas under the curve of each 
term considered separately, that is if y = f, (x) + f,(x) the area 
under this curve will be Axf,(x) + Axf,(x), since, Fig. 15, the rate 
of area production is the same whether the curves are separate or 
combined, 


VX 


x 


Figure 15 Figure 16 


G. Find the formula for the volume of a sphere. 
In Fig. 16 the sphere of radius r is cut by a plane p perpendicular 
to the x axis and moving with a velocity vx as shown. The rate of 


volume production is my*vx = x*) ux =m (r¢ux -x*vx) 


| 
c-h 
vy; = vx = y ux 
c-h | 
y = 
| 

| 

VA P 
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and the volume of % the sphere = (rx - (1/3)x9) = $,() in the in- 


ternal x =0 to x=r, V/2 = m(r?-(1/3)r?) = (2/3) 7r? and V= (¥/3)7 


vz 


h 
BJA 


S = Surface 
Q.= Volume 


Figure 17 Figure 18 


H. Find the surface of a sphere. 

If the radius of a sphere, Fig. 17, is increasing at a rate vr, 
the rate of change of volume is S ur since the entire surface is moving 
away from the center with a velocity ur. 


vQ = Sur =v. ( (4/3) = or 
=4nr? 


An alternate solution is as follows: 


In Fig. 1B the sphere of radius r and center O is cut by a plane 
moving with a uniform velocity vx, with A taken as the origin. The 
surface considered will be that produced by the circle of radius fh, at 
the intersection of the plane and the sphere, as the plane moves to 
the left in the figure. The rate of change of surface, vS = 27h vz, 


vz 
where vz is the tangential velocity. — = cse z, or vz = csc z vx 
Ux 


and h =r sin z. 
27h vz = 27r sin z*cse z vx, 
and 
sin z*cscz = wS = 2ar vz, 


and the surface generated when vx has produced a length 2r is 


2nr-2r =4nr? 


This also indicates that the surface of the portion of the sphere 
between A and the cutting plane at C is that of the convex surface of 
a cylinder of radius r and altitude AB where 


AB=r- Vp? 


{ 
Va 
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The application of the method is that of integral calculus, except 
for the derivation of the formula for the areas under the curve y = x", 
which is direct, and yields the equivalent of the integral of the 
function. 

The method of finding the equivalent of the derivatives of the 
trigonometric functions, which is unrelated to the differential 
calculus method was discovered by the writer prior to his contact 
with differential calculus. 

No attempt has been made to present this material in a form suit- 
table for presentation to high scinool students, but rather to explain 
the basic principles involved and to demonstrate the possibilities of 
the method. |Inverse trigonometric functions were not mentioned, since 
they are closely related to the direct functions. Logarithmic and 
exponential functions are not treated. No success has been attained 
in applying the velocity concept to the derivatim of the rate of 
change of these functions. 


NOTE: The integral of the sine and cosine are directly obtainable 
from Figure 10), sin x being the quantity produced by a velocity of 
cos xvx and 1 - cos x is produced by a velocity sin xvx inthe interval 
zero to x. By an extension of the method the areas under such curves 
as sin’x, sin’x, x sin x, etc., may readily be found. This method 
was developed after the above paper was submitted. 


Richland, Washington 
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MISCELLANEOUS NOTES 
Edited by 
Charles K. Robbins 


Articles intended for this Department should be sent to Charles K. Robbins, 
Department of Mathematics, Purdue University, Lafayette, Indiana. 


GRAPHICAL CONSTRUCTION TO FIND APPROXIMATE LENGTH 
OF THE CIRCUMFERANCE CF A CIRCLE 


Jerzy Cicero-Pienkowski 


Construction: 


Draw a circle. Through its center draw two straight lines (xx) and 
(yy) perpendicular to each other. Draw straight line (x,x,) parallel to 
line (xx) and tangent to the circle. Also draw straight line (y,y,) 
parallel to (yy) and tangent to the circle. This will establish three 
points (A), (B), and (J), and the square, (AOIB) whose sides are 
equal to K, the radius of the circle. From point (B) describe an arc 
with its radius equa! to four times (OA), and intersecting line (xx). 
This will locate point (C). From point (C) describe an are with its 
radius equal to twice ((A) and intersecting line (x,x,). Thus we 
shall get point ()). The distance from point (A) to point ()) will 
be approximately half of the circumferance of the circle. If we re- 
peat the same construction on the left hand side of line (yy) we will 


get point (H) and since (HA) = (AD) 


(HD) = circumferance of the circle - approximately. 
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Proof: 
By construction 
(yy) | (xx) and 
\| (yy) and (x,x,) (xx) 
(BC) = 4(0A) = 48 
(CD) = 2(0A) = 2R 


If we draw a straight line through point (C) and parallel to line 
(yy) then 


(CE) = (OA) and YX (AKC) = 90° 


Also we get two right angled triangles: (BCE) and (DCE). 
In A (BCE): 
(BE)? = (BC)?~ (cE)? 
(BE) = Vv (BC)* (cE)* = (a)? = 
= =V15R2 =k ; 


(DE)? = (pc)? ~ (ce)? 
(DE) = V(pC)? (CE)? = V(or)? — = 
= W4an?_R? = - V3. 


= 


(BD) = (BE) ~ WE) = RVI5 - = V3) 


whence 


(AD) = (AB) + (Bo) = (OA) + (BD) =R+R( VW) = 
= RO + V7I5 - V3) =R (3.14093....) 


Now the circumferance of a circle equal 27rwhere nm = 3.14159.... 


Hence one half the circumferance equal r+ 3.14159... 
Hence (4D) equal one half circumferance of a circle - approximately; 


and (HD) equal the circumferance - approximately. 


Ruffalo, New York. 
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IBN HAITHAM’ S PROBLEMS 
AND THETR GEOMETRIC SOLUTIONS 
Amir-Moez 


Ibn Haitham, (Alazen) (965-1039), has done a great deal in geometry 
and optics. He solved two problemsby constructing the roots of quadratic 
or biquadratic equations. However, he was never satisfied with his 
solutions. 

1, Problem: The circle (O) with center at O and the point P inside 
(O) are given. Find the direction of the ray through P such that after 
two reflections it passes again through P. 


Solution: Suppose PT is the desired ray (Fig. l). Clearly OT is 
the angle bisector of 7PTQ and TR, the tangent to (0) at T, is the 
angle bisector of ZM7Q. Let OP = a, OT = r, OQ= x, OR =~ y. It is 

RQ” 


easy to see r* = ay ond = —, i. . Therefore we get 


2ax? +r?x - ar? = 0. Clearly x can be constructed geometrically. 


2. The solution of Ibn Haitham’s problem: Suppose PT is the desired 
ray, (Fig. 2). Clearly 7Q, the reflected ray, is perpendicular to OP. 
If TR is the tangent to ©) at T, then T(POQR) is a harmonic pencil. 

On the other hand 7Q is the polar of R with respect to (0). Gon- 
sider the line ML, through P and perpendicular to OP. The lines AW and 
RL intersect (0) respectively at C and B. BM and CL intersect at Q. 
PQOR is a harmonic range, so is LKBR. Therefore MK passes through 0. 
Consider the harmonic pencil O(LKBR) which intersects LM at L, M, R, 
and X. Therefore LPXM is a harmonic range and this shows the con- 
struction. 

Ibn Haitham has solved the following problem by use of biquadratic. 
But the geometric construction of it is desired, 

Problem: Given A and B on the plane of the circle (0) with center 
O, find the ray through A such that its reflection on (0) passes 
through B. 


University of Idaho. 
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On A Certain Problem In Mechanics 
Murray R. Spiegel 


In a book on mechanics, a problem similar to the following appeared. 


An object having mass of 1 gram moves in a straight line on 
a horizontal frictionless plane with initial velocity of 
4 cm/sec. If it is acted upon only by a resisting force in 
dynes numerically equal to twice the instantaneous velocity 


in cm/sec, how far from the starting place does the object 
come to rest? 


It was intended that the problem be solved as follows. Letting the 
object move on the x axis starting from the origin with initial velocity 
4 cm/sec, it is clear that the differential equation for the instan- 
taneous velocity is 


(1) 


subject to the conditions 
(2) v = and 


Since (1) may be written as 


du 
(3) sane 


dx 
we have on integration 
(4) v=-2+C 
Applying conditions (2) it is seen that C = 4 so that 
(5) v=4-2 


Clearly then from (5) v = 0 when x = 2 and so the book states that the 
object comes to rest at a distance 2 cm from the starting place. 

The above problem is slightly misleading in that it is but a simple 
matter to show that the object never comes to rest. To show this, start 
with equation (1) and obtain by separation of variables and use of 
conditions (2) the velocity 


(6) v = 4e7?t 


wiich shows that v is never zero in any finite time and that as a con- 
sequence the object never comes to rest. The solution starting from 
equations (3) hides the time element in the problem and gives the false 
impression that the object does come to rest. 


Rensselaer Polytechnic Institute. 


(Will some reader explain this paradox? Fd.) 
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A Proof of the Moment-Area Theorem 
J. P. Hoyt 


In finding the deflection of statically determinate beams with 
uniform cross-sections, the following (“moment-area”) theorem is 
sometimes used: 


When a straight beam is subjected to bending, the distance 
(referred toas tangential deviation) of any point B on the elastic 
curve, measured normal to the original position of the beam, 
from a tangent drawn to the elastic curve at any other point A, 
is equal in magnitude to the product of 1/EFJ and the moment of 
the area bounded by the graph of the bending moment equation of 
the beam, the x-axis, and ordinates of the graph of the bending 
moment equation at A and B, with respect to the ordinate through A. 
(E is the modulus of elasticity of the material of which the 
beam is made and 7] is the moment of inertia of the area of the 
cross-section of the beam with respect to its neutral axis.) 


In Figure I, let arcABbe the graph of y =f (x), (where f(x), 
f'(x), and f"(x) are continuous throughout the closed interval (a,6)) 
from x = a to x = b, Let AE be tangent to AB at A and AF be parallel 
to @ . Since tangent = f'(a), EF =(b- a) fla), FD=AC=fla) 
and BD) = f(b). Hence BE = BD-FD-EF or BE = f(b) - f(a) - (b-a)f (a). 


Now although we have found the exact value of the tangential 
deviation, it 1s instructive as an illustration of summation by.in- 
tegration as well as leading to our desired result to find an integral 
whose evaluation will give BE. 

To do this, divide CD)inton equal (not necessary but sufficient) 


subdivisions of length x by the points «++ Let P 


be the point on AB whose abscissa is x ; and Q be the point on AB 
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whose abscissa is x;,;. Now x44; = x; +Ax. Let PR be tangent to 
AB at P and Qf be tangent to AB at Q. 


By the formula for tangential deviation given in the preceding 
paragraph, 


BR = f(b) fle,) - (b - f' (x) and 
BT = f(b) f(x, + Ax) - (b- x; -Ax) f (x, 
so that 
At, =TR=BR - BT 
= fla, +Qx) - -Axf' (x, + Ax) + (b-x,)[f' (x, + Ax) - f' (x,)) 
Now 


BE = At; 
for any n, so that _ 


BE = Lim At; 
t= 


Since f(x) and f' (x) are continuous throughout (a,b), we can ex- 
press f(x, + Ax) as f(x,;) + Ax + 0(A x) and f'(x, +A x) as 


f(x) + f"Cx) Ox + ofA x) where o(A x) stands for a function of Ax 
whose ratio to Ax approaches 0 asAx approaches 0. 


Hence upon substituting and simplifying, 


At, = (b- x) o( Ax) 
And 


BE = lim At, = lim [(b - x) f" @)A «) + Lim 
n-@ 
1= 


But 


o( Ax) = 0. 


So that t= 


b 
BE = lim (b-x,) f"(x,) Ax = f (b - x) f"(x) dx 
ts ° 


by the fundamental theorem of the integral calculus. 
This last integral can be recognized as the moment of the area 


bounded by y = f(x), y =0, x =a, x = 6, with respect to the line 


o( Ax) 
i= 
Lim 
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x = 6. Since the second derivative with respect to x of y inthe ea- 
uation of the elastic curve y = g(x) is 1/E] times V,where Mis the 


the bending moment equation of the beam, the moment-area 


theorem 
follows. 


Apart from the usefulness of this theoreninthe strength of materials, 
it merits inclusion in calculus texts (I have never seen it in one) for 
two reasons: one, in deriving it we use the integration process to 
find the limit of a sum whose limit we already know: two, it can also 
be used to find not-so-easily calculuable moments from easily cal- 
culable tangential deviations. Such problems are also useful to an 
understanding of the theorem. One such example would be: Find the 
moment of the area lying in the first quadrant bounded by the lines 
y = sin x, y= 0, x= 7/2. This moment would be given by the integral 


m/2 
f (n/2 -x) sin xdx = 7-1, 
‘ 2 


This involves integration by parts. To apply the moment area-theorem, 
find the simplest function, -sin x, whose second derivative with re- 
spect to x is sin x. Draw the tangent to the graph of y =~ sin x at 
x = 0, and find the tangential deviation at x = 7/2. This gives the 
required moment. From figure II, the steps are obvious and simple. 


tan ZDOA = -cos 0 


AD = —0OA 


CD =AC ~AD= 


Figure II 


. Naval Academy. 
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CURRENT PAPERS AND BOOKS 


Edited by 
H. V. Craig 


This department will preere comments on papers prarienedy published in the 
MATHEMATICS MAGAZINE, lists of new books, and book reviews. ‘ 

In order that errors may be corrected, results extended, and interesting 
emporte Meshes illuminated, comments on published papers in all departments 
are invited. 


Communications intended for this department should be sent in duplicate 


to H. V. Craig, Department of Applied Mathematics, University of Texas, 
Austin 12, Texas. 


Principles and Techniques of Applied Mathematics. By Bernard 
Friedman. John Wiley & Sons, 440 Fourth Ave. New York 16, N.Y. 315 
pages, $8.00. 


A unique combination of abstract theory and practical application 
marks this book. It is the latest addition to Wiley’s Applied Mathe- 
matics Series, edited by I.S. Sokolnikoff. 

Deploring the widening gap between pure and applied mathematics, 
Dr. Friedman states his belief that this gap is an illusion. He con- 
tends that the study of abstract systems helps solve concrete problems, 
and that the studyof specific problems may, in turn, suggest interesting 
generalizations to the pure mathematician. Principles and Techniques 
of Applied Mathematics therefore sets out to show how the powerful 
methods developed by the abstract studies can be used to systematize 
the methods and techniques for solving problems in applied mathematics. 

Two main themes carry out Dr. Friedman’s approach. The first is a 
demonstration of how the abstract theory of linear operators can be 
used to unify the techniques of applied mathematics. Secondly, Dr. 
Friedman develops and explains specific techniques that can be used to 
obtain explicit solutions of partial differential equations. 

Among the novel features appearing in this volume are a presen- 
tation of the Schwartz theory of distributions and the S-function 
with applications; a simple proof for the reduction of a matrix to 
Jordan canonical form, with the proof using generalized eigenvectors; a 
simple technique for obtaining the complete eigenfunction expansion 
for a second-order ordinary differential operator; a formal method for 
solving partial differential equations which includes all transform 
methods such as the Laplace and Fourier transforms; and the useof the 
general theory of the linear operators to obtain the Fredholm theory 
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of integral equations. Admitting that his emphasis is on the techni- 
ques, Dr. Friedman maintains that a logical analysis of a problem 
leads automatically to the proper methods for its solution. 


Richard Cook 


Symposium on Monte Carlo Methods. Edited by H. A, Meyer. John 
Wiley & Sons, 440 4th Ave., New York 16, NY. 382 pages, $7.50, 


The first full-lenth treatment of this subject, the new book is 
the work of twenty-two leaders in this promising field, writing about 
their own research and applications. 

The volume is the outgrowth of the symposium conducted by the 
Statistical Laboratory of the University of Florida, and sponsored 
by the Wright Air Development Center. It is the latest addition to the 
Wiley Publications in Statistics, Walter A. Shewhart and ~~ S. S. Wilks 
Editors. 

“The possibility of going from a physical process directly to a 
numerical solution without benefit of an intermediary differential 
equation,’ Dr. Meyer states, “is such an unorthodox and revolutionary 
approach that the imagination is kindled to reconsideration of un- 
solved problems of old.” 

Originally refined for solving atomic energy problems, Monte Carlo 
methods have begun to spread into other specialties such as market 
analysis, quality control, public-opinion polling, and genetics. 

The papers appearing in this book range from the simplest to the 
highly theoretical and include the following: An Introductory Note by 
A.W. Marshall; Generat.wn of Pseudo Random Numbers by Olga Taussky 
and John Todd; Phase Shifts - Middle Squares — Wave Equations by 
N. Metropolis; A General Theory of Stochastic Estimates of the 
Neumann Series for the Solution of Certain Fredholm Integral Equations 
by G.E. Albert; Neighbor Sets for Random Walks and [Difference Equa- 
tions by Theodore S. Motzkin; Monte Carlo Computations by Nancy M. 
Dismuke; Applications of Monte Carlo Methods to Tactical Games by 
S. Ulam; Conditional Monte Carlo for Normal Samples by Hale F. Trotter 
and John W. Tukey; Monte Carlo Techniques in a Camplex Problem about 
Normal Samples by Harvey J. Amold, Bradley D. Bucher, Hale F. Trotter, 
and John W. Tukey; and An Application o the Monte Carlo Method to a 
Problem in Gamma Ray Diffusion by Martin J. Berger. 

The balance of the papers include: Stochastic Calculations of Gamma 
Ray Diffusion by L.A. Beach and R.B. Theus; Application of Multiple 
Stage Sampling Procedures to Monte Carlo Problems by A. W. Marshall; 
Questionable Usefulness of Variance for Measuring estimate accuracy 
in Monte Carlo Importance Sampling Problems by John E. Walsh; Expert- 
mental Determination of Eigenvalues and Dynamic Influence Cooefficrents 
for Complex Structures such as Airplanes by C. W. Vickery; Use of 
Different Monte Carlo Sampling Techniques by Herman Kahn; A Theoretical 
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Comparison of the Efficiencies of Two Classical Methods and a Monte 
Carlo Method for Computing One Component of the Solution of a Set of 
Linear Algebraic E;aations by J. H. Curtiss; A Description of the 
Generation and Testing of a Set of Random Normal Deviates by E. J. 
Lytle, Jr.; Machine Sampling from Given Probability Distributions by 
James W. Butler; A Monte Carlo Technique for Obtaining Tests and Con- 
fidence Intervals for Insurance Mortality Rates by John E. Walsh; and 
Experiments and Models for the Monte Carlo Method by Alwin Walther. 

In addition, an extensive bibliography gives abstracts and sum- 
maries of existing | .«rature, many of them formulated in the authors’ 
own words. A consicerable portion of this literature was previously 
available only in the form of government technical reports which were 
frequently of limited distribution. The appendix also contains a sec- 
tion on random digits, a section dealing with empirical sampling, and 
a selection of papers making use of stochastic processes, including 
a bibliography prepared by Dr. David G. Kendall. 


Richard Cook 


Irrational Numbers. By Ivan Niven. John Wiley & Sons, New York 
164 pages. $3.00. 


Number Eleven of the Carus Mathematical Monographs, /rrational 
Numbers was published in August by the Mathematical Association of 
America. It will be distributed by John Wiley & Sons. The author is 
Ivan Niven, professor of mathematics at the University of Oregon. 

This book provides an exposition of some central results on irra- 
tional, transcendental, and normal numbers. There is a complete treat- 
ment by elementary methods of the irrationality of the exponential, 
logarithmic, and trigonometric functions with rational arguments. 

The approximation of irrational numbers by rationals, up to such 
results as the best possible approximation of Hurwitz, is also given 
with elementary techniques. The last third of the monograph treats 
normal and transcendental numbers, including the transcendence of 7 
and its generalization in the Lindemann theorem, and the Gel fond- 


Schneider theorem. 
Richard Cook 


Digital Differential Analyzers. By George F. Forbes, 10117 Bartee 
Avenue, Pacoima, California, Paper cover, approximately 200 pages 
including appendices and nine page index. $7.50 


An applications manual dealing with mathematics in the form speci- 
fically related to digital differential analyzers. A knowledge of 
differential and integral calculus is presupposed. Familiarity with 
more advanced mathematics and with a specific DDA is desirable but 
not essential. 
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The common functions are discussed in some detail, as well as servo 
techniques and limiting methods. Differentiation, the more advanced 
functions, simultaneous equations, and implicit algebraic functions 
are described in a more general way. Following a section on errors of 
the DDA, the more advanced sections of the text deal with complex 
functions, coordinate system transformations, conformal mapping, 
partial differential equations, and curve analysis. 

Several illustrative problems are discussed. The MDA is treated 
both as a mathematical device and as an engineering and scientific 
tool. 


This year International Business Machines estimates that 
7,500 Mathematicians will be needed to man Computers now on 
order, of which about 1,500 should be Ph. D’s. About 250 
will graduate this year! 


From Aviation Week, Sept. 17, 1956, p. 25 
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PROBLEMS AND QUESTIONS 
Edited by 
Robert E. Horton, Los Angeles City Gllege 


Headers of this department are invited to sumbit for solution problems 
believed to be new and subject matter questions that may arise in study, 
in research, or in extra-academic situations. Proposals should be accom- 
panied by solutions, when available, and by any information that will assist 


the editor. Ordinarily, problems in well-known textbooks should not be 
submitted. 


Solutions should be submitted on separate, signed sheets. Figures should 
be drawn in India ink the size desired for reproduction. 

Send all communications for this department to Robert £. Horton, Los 
Angeles City College, 855 N. Vermont Ave., Los Angeles 29, California. 


PROPOSALS 
285. Proposed by Maj. H.S. Subba Kao, New Delhi, India. 


It is well known that given any two areas enclosed by closed curves 
in a plane there exists a straight line which bisects both the areas. 
Is it possible to construct such a line when the two areas are: a) a 
triangle and a parallelogram, b) two triangles? 


286. Proposed by Joseph Andruskiw, Seton Hall University. 


Evaluate 


n/2 
l+sin“x 


287. Proposed by Chih-yt Wang, University of Minnesota. 
Show that 
cot(7/11) - cot( 2/11) + cot( 47/11) - cot(87/11) + cot( l6n/11) =v 11 


288. Proposed by Norman Anning, Alhambra, California. 

A particle moves on a Fermat Spiral, r? = a6, under a central 
force directed to the origin. Show that the law of force involves the 
seventh power of the radius in a denominator. 


289. Proposed by J.M. Howell, Los Angeles City College. 

Prove or disprove the following statement: The diagonal of a rec- 
tangular parallelopiped can be any odd integer other than one or 
five, and the edges relatively prime integers. 

290. Proposed by Leon Bankoff, Los Angeles, California, 


The hypotenuse AC of a right triangle ABC is divided into unequal 
segments x and y by the point of contact of the incircle. Given seg- 
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ments x and y: a) construct the triangle ABC and its incircle, b) with- 


out using the fourth proportional, upon x as a base, constructa rec- 
tangle equivalent to triangle ABC. 


291. Proposed by M.S.Klamkin, Polytechnic Institute of Brooklyn. 
Determine the minimum of 


where x, > 0 


Solutions 


Probability of Convexity of Four Points 


264. { March 1956] Proposed by Norman Anning, Alhambra, California 


Four points are thrown at random on a plane. What is the probability 
that they will be the vertices of a convex quadrilateral? 


Solution by Huseyin Demir, Kandilli, Bolgesi, Turkey. The probability 
of three points to form a triangle is evidentlyl. Let then ABC be the 
triangle formed by three of the points. For the points A, B, C, |) to 
form a convex polygon it is necessary that the fourth point () be in 
the region of ABC where we have excircles. These regions of the plane 
correspond to the angles 4, B, C of the triangle ABC. Hence 


p = (A +B + C)/whole plane = 7/27 = 


Also solved by Richard K. Guy, University of Malaya, Singapore; 
John M. Howell, Los Angeles City College and the proposer. 


A Volume of Revolution 


265. [ “arch 1956] Proposed by Stephen Armstrong, Union College 
Schenectady, New York. 


What is the volume of revolution of the ellipse 


revolved about the line x = y? 


I. Solution by Lawrence A. Ringenberg, Eastern Illinois State 
College. Assume a > b > 0. Let S denote the sector with central angle 
7/2 which is bounded on the left by the lines y =+x and on the right 
by the given ellipse. If (x,y) is a point in S, its distance from the 


>> 
é 
x y 
! 
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line y = x is (x -y)/ 2. Since the desired volume V is twice the 


volume generated when S revolves about the line y = x, we may compute 
it as twice the double integral of V/27(x-y) over S. Thus 


V =2V2 ff(x- y) as 
Set 
ab 
+b? ‘ 


A fly) 
(x-y) dx dy+2 /27™ f (x-y) dx dy 
y 


4V2 
3Va2 + b2 


Comment. If the ellipse is rotated through one-half a revolution, the 
volume V,. generated may be computed as the double integral of 
¥27 (x -y) taken over that portion of the xy-plane which lies within 
the ellipse and below the line y = x. It may be interesting to note 
that V, could also be computed as follows. Let the lines y =+x divide 
the interior of the ellipse into four quadrants S,, S,, S,, Sy 
numbered counterclockwise starting with S, = S, S of the preceding 


paragraph. Let V;, « = 1,2,3,4, be the volume generated when S, rotates 
one-half revolution about y = x. Then if V is the volume obtained in 
the foregoing solution and 


4 V2 nab? 
, Obtained from V by interchanging a and 6, we have 


3 Va? + b? 
V, = V/4, V,=V¥,=V/4, and Vo = 


II. Solution by C.N, Mills, Augustana College, Sioux Falls, South 
Dakota. There are two equal generating surfaces revolving about the 
axis y =x. ‘fhe area of each surface is 


ab [ b ] 

— - 2 arc sin 

: Va2 + b? 

The distance of the centroid of each surface from the axis is 
2 a? b 


3 + b? 


Applying the theorem that the volume of a solid of revolution is equal 
to the area of the generating surface times the distance through which 
the centroid moves, we obtain 


[- - 2 are sin 


0 
v= 2v2 of f 
-A 
and 
| 
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4 
| 
3 Va? + 


Also solved by Huseyin Demir, Kandilli Bolgesi, Turkey; J.M.C. 
Hamilton, Los Angeles City College; Richard k. Guy, University of 
Malaya, Singapore; Sister M. Stephanie, Georgian Court College, New 
Jersey and the Proposer. 


Points Inverse in a Circumcircle 
266.(March 19956] Proposed by Huseyin Demir, Zonguldak, Turkey. 


If M and M’ are points inverse to each other with respect to the 
circumcircle of a triangle ABC then prove that: 


2 BMC + ZBM'C = 2ZA 
Z CMA + ZCM'A = 22B 


2 AMB + ZAM'B =22C 


I. Solutio by Richard K. Guy, University of Malaya, Singapore. 
In triangles OOM and ‘OC angle O is common and as OV + OM = OC? we 


have a = on” Hence the triangles are similar and {QW'C = JOM, 
In the same way JOM'B = ¥ (BM. Adding these to YOWC and JOMB we 
have YBMC +ZBM'C = - {COM + - YBOM = 2 YA. 


Simmlarly we have {CMA +YCM'A = 24B and JAMB = JAM'B = 


II. Solution by Maimouna Edy, Hull, PQ, Canada. Represent points 
A, B, C, M, M’ by complex numbers 2, z' respectively. 
Let parentheses represent cross ratios and the bars the complex con - 
jugate. We then have: 


() (25, z) = (zy, By z’) 


This says that the homographic transformation which sends z,, z, 2 
into 1, 0, © respectively, that is the transformation of the given 
circle into the axis of reals, sends z andz’ into two conjugate complex 
points. 

Now equation (1) reads explicitly 


Evidently 
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Therefore 


z,-2 272, 
This means that for oriented angles, 
(angle MC, MB) + (angle MC, M’B) = 2(angle AC, AB). 


The oriented angles form an additive group isomorphic with the multi- 
plicative group of the unit circle. In other words, we may take 
arbitrary measures of our angles and add them (mod 27). The other two 
relations are proven similarly. 

Bankoff’s solution also noted the necessity for proper orientation 
of the angles. 


Also solved by Leon Bankoff, Los Angeles, California; J.W. Clawson, 
Collegeville, Pennsylvania; and the pro pser. 


Liwit of a Radical Sus 


267. (March 1956] Proposed by Alan Wayne, Cooper Union School of 
Engineering, New York. 


n 
If (k -1) , prove that lum 


k=1 


I. Solution by Paul Schillo, University of Buffalo. Let P(c,u,r) 
denote any polynomial in u of degree t, where 0 < t <r and where 
c is the coefficient of u’. If p is a positive integer and 6 7 0, then 
lettingQ.(1/m) denote a function f of m and k such that mf(m,k) is 
bounded for all positive integer pairs (m,k), in which m is greater 
than the greatest integral root of P(b,m,r + p) = 0 and k < m,it 
follows (by the fundamental theorem of the integral calculus) that 


a 1/p 1/p 


poy P(b,m,r + p) 


a 1/ Pp a 1/p 
= "| = ~) 
b b 


In the given problem, k” - (k - I)" =P(n, k, n - 1) and a"**= 


Pl, a, n+ D, so Lim A, 
It might be added that the above limit holds if the polynomials have 
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complex coefficients, since it is the line integral of an everywhere 
analytic function. 


II. Solution by the proposer. Make use of the (u, v) rectangular 
plane coordinate system in which u = x + iy, v =x - iy, andi = V-L 
In this (u, v) plane, the distance between, the two ppints P (u,v, ) 
and Pi(u,, v,) is given by P,P, = (u, - u,J*(v, - v,)?. Moreover, the 
length of the arc of any rectifiable curve v = f(u) joining P, and 


P, is given by ‘, 
PP, (dv /du)” du. 


u 
1 
Now consider the curve v = u", with its points 
Py(k/m, (k= 0, 1, 


Clearly for k>0, I" / Then 


2 


Similar limits of sequences of series can be evaluated by this 
“transformation” method. 

Also solved by Maimouna Edy, Hull, P.Q. Canada and Chih-yit Wang, 
University of Minnesota. 


Coaxal Circles 


268. (March 19956] Propsed by J.W. Clawson, Collegeville, Pennsylvania. 


Three coaxal circles, centers at A, B and P have the common points 
C and PD. Any straight line is drawn through C cutting the circles 
again in the points L, M and N respectively. Prove that the ratio 
IN/NM equals the ratio AP/PB. 


Solution by Chih-yi Wang, University of Minnesota. Let the feet 
of perpendiculars from A, B and P onto the line LCMN be A’, B’ and P’ 
respectively. Then we have 


CP’ = CN/2, CB' = CM/2, A'C = LC/2; 


PB. BP AB AB+ BP AP 


or 


BS6) 
Hence 
1 
lim = PP (nu®* du = 2VnAn+ 
0 
P'B' PB’ 
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whence LN _ EN A‘C AP 
NM NM/2 CB’ + P'C’ PB 


Remark: If the line passing through C is the line CD, L, M and N 
coincide and in this exceptional case the ratio LN/MN is of the form 
0/0. 

Also solved by Huseyin Demir, Kandilli, Bolgesi, Turkey; Maimouna 
Edy, Hull, P.Q.,Canada; Richard K. Guy, University of Malaya, Singapére; 
and the proposer. 


The Product of Power Series 


269. March 1956] Proposed by M.S. Klamkin, Polytechnic Institute of 
Brooklyn. © 


,a-! Ree 
Find the sum >> + 
3! 


Solution by Maimouna Edy, Hull, P.Q., Canada. 


@ 
n 


This result is valid for |x| < L 


Also solved by L. Carlitz, Duke University; V.D. Gokhale, Atlanta 
University, Georgia; John M, Howell, Los Angeles City College; Alice 
A. Huffman, University of California at Rivers ide ;Lawrence A. Ringenberg, 
Eastern Illinois State College; Paul Schillo, University of Buffalo; 
Chth-yt Wang, University of Minnesota and the proposer. 


Circles In A Crescent 
270. (March 1956] Proposed by Leon Bankoff, Los Angeles, California. 


A maximum circle is inscribed in a crescent formed by a semicircle 
and a quadrant of a circle. Find a general expression for the radii 
of consecutively tangent circles touching the sides of the crescent, 


| 
n! 
| 
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the first touching the maximum circle, the second touching the first 
and so on. 

Solution by Huseyin Demir, Kandilli, Bolgesi, Turkey. Let the 
given circles (0), (0') intersect each other at A and B, and let the 
center and radius of the nth circle be denoted (0), 7, respectively, 

We invert the figure with center at A, k* = AB? = 4R? being the 
power. Under the inversion, (0) is inverted into its tangent line 


BH,, and (0') into the line BO’, forming an angle of 2a= 45°. The 
circles. (,), inverse of (0;), from a series of tangent circles in- 


scribed in the above angle. Let ({)2,) touch Bly at H,. Then we may 
easily find that the radius p,=1,H, of (2 ,) is given by 
1- sina 


p = 
1+ sina 


where 
Po = Ndlo = BH, tan a = 2R tan a, 


Drawing the common tangent AT,7T) to the inverse circles (0), 
(Q,) we write 


AT, = AT, /ATA? = /(AQ? ni?) 
Denoting the projection of {, on AB by K, we have 
2 2- 
4, =m, * 
(2k + p,)? + BH? - p2 
4h? + 4R p, + (p, cot a)? 


tana 


AR? + 4R + pe cot? a 


Substituting the value of P, in the above expression we arrive at the 
desired result, namely 


kK 
1+ + sin +(1 - sin 7/8) 2") sec*" 8 cot 1/8 


Also solved by J.W. Clawson, Collegeville, Pennsylvanta and the 
proposer. 
TRICKIES 
A trickie is a problem whose solution depends upon the perceptian of the 


key word, phrase or idea rather than upon a mathematical routine. Send us 
your favorite trickies. 
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T 24. A certain physical society is planning a ballot for the election 
of three officers. There being 3, 4, and 5, candidates for the three 
offices, respectively. There is a rule in effect (in order to eliminate 
the ordering of the candidates on the ballot as a possible influence 
on the election) that for each office, each candidate must appear in 
each position the same number of times as any other candidate. What 


is the smallest number of different ballots necessary? [Submitted by 
M.S. Klamkin} 


T 25. Given the sequence 1, 1, 2, 3, 5, 8, ... where @; =a,.,°+ a, , 
find an expression for the nth term, a,, not in terms of a recursion 
formula, Generalize the sequence so that a, and a, may be any two 


numbers. [Submitted by Ben B. Bowen), 
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QUICKIES 


From time to time this department will publish problems which may be 
solved by laborious methods, i which with the proper insight may be dis- 
posed of with dispatch. Readers are urged to submit their favorite problems 
of this type, together with the elegant solution and source, if known, 


Q 183. What are the ratios of the volumes and surfaces of the inner 
“half” to the outer “half” of the torus generated by revolving the 
x24 (y - [Submitted by V.C. Harris} 


@ 184. Simplify 


sin x + sin 2x+ sin 3x +... + Sinnz 


cos x +cos 2x+ cos +... + cos nx 
[Submitted by M.S. Klamkin] 
Q 185. Prove that 
ner 
where s <k - ] 
r=0 
[Submitted by M.S. 
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PERCENT TABLE 


G.0.JAMES ENTERPRISES 
5209 Windermere Ave. 
Les Angeles, 41, Calif. 


ANNOUNCES 


THREE % TABLES DESIGNED FOR SAMPLING 


In penphlet form, gives percentages for numerators to 
50, (150 with only one additional step), and any de- 
nominetet ~ percentages accurate te .1% from 0 ¢04.8% 
os very goed approximation from thet point on. 


Quick epproxiaate % Teble, in menile folder size, covers 
to and any denominator. 


Quick % Table, in manile folder site, Covers numerators 
to 3, denomtaetors in 10's or 25's to 830, percentages 
to .&%, 


ALL THREE TABLES FOR OWLY 
$3.50 


TWICE AS FAST AS 
A SLIDE RULE YET 
cneaper ! 
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